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ABSTRACT 


The primary analyses performed on data collected during low- 
level tracer experiments conducted over metropolitan St. Louis, Mis- 
souri, are described. Values of dispersion parameters derived from 
the tracer data are related to readily derived or measured meteoro- 
logical indices of turbulence. The results are graphically presented 
in terms of best-fit curves as functions of downwind distance and 
travel time. Comparisons are made with the results of previous dif- 
fusion experiments conducted over relatively uncomplicated terrain 
in open country. It is concluded that for low-level point sources the 
urban area affects cross-wind dispersion primarily by enhancing the 
initial size (i.e., close to the source) of the plume. As the plume 
becomes much larger than the size of eddies created by the local ob- 
structions, the dispersion approaches that associated with flow over 
open country. In the vertical, significantly enhanced dispersion as 
well as an enlarged initial spread of the tracer cloud seem to occur; 
this enhancement in the rate of dispersion over that in open country 
appears somewhat greater for stable than for unstable meteorological 
conditions. Values of dispersion parameters from limited tracer ex- 
periments in other urban areas are similar to those reported here for 
St. Louis under the same gross meteorological conditions. 
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ST. LOUIS DISPERSION STUDY 
VOLUME П- 
ANALYSIS 


INTRODUCTION 


The St. Louis Dispersion Study consisted of a series of experi- 
ments in which fluorescent zinc cadmium sulfide particles were releas- 
ed and traced across the relatively flat urban area of metropolitan St. 
Louis. The measurement program began in the spring of 1963 and 
ended in the spring of 1965. Over this 2-year period 26 daytime and 
16 evening experiments were conducted in seven series. Dissemina- 
tion of fluorescent particles from either of the two pre-selected sites 
near ground level was usually 1 hour long. For each experiment, 
measurements of total dosage at the surface were obtained on three 
nearly circular arcs at distances between 1/2 and 10 miles from the 
dissemination site. The use of particular arcs varied, usually depend- 
ing on appropriate forecasts of wind direction and speed. In addition, 
time-sequential measurements of dosage at the surface were made on 
each arc at a few locations near the anticipated mean centerline of the 
tracer cloud. During nine experiments, total dosage was measured at 
several heights along the tether of a balloon flown at a single site. The 
site was usually a park or vacant lot between the inner and quter sam- 
pling arcs and as near as possible to the anticipated mean centerline. 
Also, a mesometeorological network composed of three stations on the 
periphery of the urban area and an instrumented television tower 
(КМОХ) in the downtown area provided continuous records of wind, 
temperature, and relative humidity. The TV tower was instrumented 
at three levels to provide information on the vertical gradients of wind 
and temperature. Single-theodolite (pibal) observations of winds aloft 
and measurements of winds near the surface were made at the tracer 
dissemination site; free or tethered radiosonde ascents were made from 
the roof of a building in the downtown area; and transponder -equipped 
tetroons were released near the dissemination site. The tetroons were 
tracked by radar located at Lambert Field, northwest of the metropoli- 
tan area. 


Volume I of this report (APTD-68-12) presents a detailed descrip- 
tion of the experimental equipment and procedures employed in the 
study, and provides in tabular form the dispersion and related meteoro- 
logical data collected during the tracer experiments. Because Volume 
I is mainly data tables, it is not being given general distribution. 


Volume II describes the primary analyses performed on the dis- 
persion and related meteorological data collected during the tracer 
experiments. The main objectives of the analyses were to obtain at 
least gross estimates of (1) the values of dispersion parameters over 
urban areas and (2) the relation of these parameters to meteorological 
indices of turbulence. A secondary objective was to compare the re- 
sults of the St. Louis Dispersion Study with those of past experimental 
programs, especially those conducted over relatively ''ореп'' country. 
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DESCRIPTION OF TRACER DATA 


TOTAL DOSAGE AT SURFACE 


For our purposes, measurements of total dosage were converted 
to equivalent concentration, X (particles/m"), by the equation 


X SEa „5 
(v) (e) (ФӘ 
where: 
n = total number of particles (dosage) 
v = flow rate of sampler 
e = collection efficiency of sampler 
Dt = time duration of tracer dissemination 


For each sampling arc of an experiment, values of X were plotted as a 
function of azimuth direction (degrees) from the tracer dissemination 
site. A continuous, smoothed curve was drawn through the plotted 
values. Whenever necessary, curves were extrapolated to zero соп- 
centration. In such extrapolations, the continuity and similarity of 
distributions between sampling arcs were maintained, and, whenever 
feasible, curves or portions of curves were made to resemble a Gaus- 
Sian shape. 


Multi-peaked or otherwise complex cross-wind distributions of 
tracer material occasionally occurred, particularly in data from the 
close-in sampling arcs. Channeling of the airflow due to the spatial 
distribution of obstructions or locally induced circulations may be 
largely responsible for this physical appearance of the distributions. 
Effects of sampler exposure were considered to be insignificant, since 
the placement of samplers at the nominally designated locations was 
somewhat random over the 2-year period of the experimental program. 


SEQUENTIAL DOSAGE AT SURFACE 


Measurements of sequential dosage at the surface were trans- 
formed into concentration X, by the formula 


n 


"s ^ (v) (е) (Dj 


where: 
number of particles (dosage) 
flow rate of sampler 


4 8 
Ho 


е 
D: 
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collection efficiency of sampler 
time duration of sequential sampler interval 


For each designated location, a plot was made of X, as a function of 
time. The physical appearances of the resulting histograms differed 
somewhat for daytime and evening experiments. 


For daytime experiments, the histograms showed that the number 
of major concentration peaks decreased with downwind distance from 
the tracer dissemination site. At inner sampling arcs, concentration 
was usually continuous with time, though varying widely. At interme- 
diate arcs, the sequential concentration often went to zero, suggesting 
the occurrence of ''puffs'' of tracer material. At greater downwind 
distances, the concentration patterns usually became more uniform in 
appearance. From these patterns, it is inferred that the horizontally 
meandering cloud is carried selectively aloft by convective motions. 
When most of the cloud has been affected by these motions, convective 
overturning results ina relatively uniform vertical distribution of 
tracer material. 


For evening experiments, the histogram patterns generally were 
more uniform, although many suggested the existence of weak convec- 
tive activity. Temperature gradient data from the TV tower, analyzed 
by the authors and by Schiermeier (1967), also suggest the occurrence 
of such activity, especially in the lowest 250 feet above the ground. 
The vertical extent of convective activity probably depends strongly on 
details of the vertical temperature structure since air flowing into 
urban areas from the surrounding countryside must always be ina 
transitional state. 


TOTAL DOSAGE IN THE VERTICAL 


During nine of the experiments, total dosage was measured at 
several heights, to a maximum of roughly 1000 feet above ground, 
along the tether of a balloon flown аға single location. Passage оға 
squall line during one of these experiments and of a front during an- 
other experiment prevented these measurements from being used to 
develop generalizations. Of the remaining seven experiments (four 
daytime and three evening), the measurements of one daytime and one 
evening experiment showed an anomalous increase of dosage with 
height above the surface. А possible explanation of the anomalies 15 
that on both occasions considerable vertical shear of wind direction 
occurred in the lower atmosphere. The sampling sites were far from 
the mean centerlines of the tracer clouds, but with respect to the dis- 
semination site were in the azimuth direction of this shear. The 
remaining measurements were too limited to allow direct computations 
of parameter values or to provide definitive information concerning 
concentration distributions. 
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DISPERSION PARAMETERS 


For the purposes of this report, 32 of the 42 experiments yielded 
usable dispersion data; 22 of these were conducted in the daytime and 
10 in the evening. Many of the excluded experiments yielded usable 
data for which analysis could not be made in the conventional manner 
used here. A separate publication is planned to report analyses of 
these data. 


NORMALIZED AXIAL CONCENTRATION 


From the plots of equivalent concentration at the surface, X, 
(particles /meter?) versus azimuth (degrees), 0, from the tracer dis- 
semination site for each sampling arc of an experiment, peak or axial 
values, X,, were determined. These values were converted to nor- 
malized peak concentrations, (X/Q), (seconds/meter3), through divi- 
sion by the tracer emission rate, Q (particles/second). 


CROSS-WIND PARAMETERS 


From each plot of X versus 0, values of X were read at equally 
spaced azimuth intervals, beginning at an edge of the tracer cloud. 
The intervals were equivalent to the tracer measurement resolution. 
Measures of the statistical cross-wind standard deviation, 9, (meters), 
and the cross-wind integrated concentration, CIC (particles /meter?), 
were then computed from these values. 


The equation for the root-mean-square standard deviation, 7, 
rewritten from Panofsky and Brier (1963), was modified to allow cal- 


culations of 7.;: 1/2 
y г 
М Zf.(L) - ( ХІ, 1) 
_ i i i i A 
с = (1) 
N2 

where: 

4 = frequency class interval 

I; = deviation of class midpoint i from the assumed mean 


in terms of class intervals 
f; = frequency of the distribution for class interval i 
N = total fj, i.e., 24; 


When а plot of X versus 0 is considered to Бе а frequency histogram of 
X in terms of 0, the values of X at the equal azimuth increments may 
be taken as the mean values of increments for which they are the mid- 
points. Thus, Xi = fi and ХХ: = М. The proper A is not Ag (azi- 
muth degrees) but rather ây (arc length) written 


гу = (Ag) (X) («/180) (2) 


XX.X. 
where X, mean arc radius, is defined by—}— and X; is the equivalent 


хх. 


concentration at a given sampling site along the arc located а distance 
Xj from the dissemination site. 


With these definitions equation (1) may be rewritten as 


1/2 


2 
хх. (x. I. ) - (zx. 1.) > 
ту = > Ay e -0- (3) 
(>х,) j 


Surface cross-wind integrated concentration, СІС,, defined as 


s? 


f X(y)dy was approximated by ХХ; А When 8, from equation (2) is 


у: 
substituted into the latter, the appropriate formulation becomes 


УХ.хХ. 
сіс, =) ( X) (49) с ) (+35) | (4) 


VERTICAL PARAMETERS 


Estimates of vertical dispersion parameters were based on the 
surface tracer measurements and hence represent only "effective" val- 
ues. As previously noted, the measurements of tracer material in the 
vertical were too limited for direct computation of vertical parameter 
values. The rationale for describing dispersion in the vertical differed 
for daytime and evening experiments. 


Daytime Experiments 


For each arc of a daytime experiment, a derived estimate of the 
vertical standard deviation, oz (meters), was computed from the appro- 
priate mass continuity equation, assuming no loss of tracer in transit, 
uncorrelated horizontal and vertical tracer distributions, and a Gaus- 
sian distribution of material in the vertical. This equation may be 


written as 
Q= | | а x(y Х(2) dy d (5) 
f. dibus 


quantity of tracer material disseminated per unit time 
- effective transport wind speed 


where: 


а! O 
| 


Since u is assumed to be independent of y and z and foy is CIC,, 


Q _ со 
acc C / чаа 
5 о 
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With Х(2) = 1/2 exp (-22/222), i.e., the function for a Gaussian dis- 
tribution with perfect reflection of the tracer material at the ground, 
the appropriate integration yields 


Q Vin 
woe “2 7 (6) 


E 72 = (CIC,) (а) 4/2”) 


Estimates of the effective transport wind speed, U, were obtained 
from mean tetroon trajectories, except when the surface projected posi- 
tion of the tetroon differed significantly from the centerline of the tracer 
cloud or when the elevation of the tetroon was beyond the top of the mix- 
ing layer that was assumed to contain the bulk of tracer material. On 
these other occasions, values of ü from mean wind profiles measured 
at the TV tower or by pilot balloon ascents were substituted. In experi- 
ments for which suitable sequential dosage samples were obtained, the 
measurements of u based on observed wind or tetroon data, especially 
the latter, agreed well with t values based on the timing of tracer cloud 
passages across sampling arcs. Estimates of travel time, t, of the 
tracer cloud to various arcs were made from 


t-x/u (7) 


where: 
X - mean arc distance from dissemination site 


Evening Experiments 


For the evening experiments, which were conducted under more 
stable meteorological conditions, the vertical profiles of wind and tracer 
material were assumed to be significantly interrelated. In view of this 
assumption and of the limitations of the tracer data, a simple box 
model was derived to describe dispersion in the vertical. Basically, 
the model assumes that the tracer material is uniformly mixed to some 
height h. The depth of the uniformly mixed layer, h, was assumed to 
increase linearly with travel time, t, of the tracer cloud between arcs. 
That is, 


h, = h, T K (ts - t.) (8) 


2 1 
where the subscripts refer to the particular arcs and К |2 is a constant 
for these arcs for a particular experiment. The wind profile was as- 
sumed to be of the form uz; = az where z is height and а and n are 
constants for each experiment. These particular assumptions appeared 
reasonably to characterize the appropriate tracer and meteorological 
measurements. 


For the simple box model, equation (5) reduces to 
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О = (CIC) (a, ) (h) (9) 


or h = Q/CIC. В, 


where йр is the effective mean transport wind through depth, h, аға 
particular arc. The appropriate uj for h at each arc was defined by 


the relation E 


| X(z)u, dz 
Бо 


Чч ---- (10) 


| X(z) dz 
о 


where X(z) in this case is a constant, С. 


Integration of this equation with the substitutions ч, = az? and 
Х(2) = С yields 
n 
ah 
ч, = ——— 11 
"h п + 1 Of) 


When this value for üj is substituted into equation (8), the relation 
becomes 


1 


[n + 1 Q n + 1 
Таға n 


Now it is necessary to determine horizontally and vertically 
averaged travel speeds, Up, based on the horizontal variation of verti- 
cally averaged speeds, uy. Since h from equation (8) has been defined 


as a function of t, 
2 
| а (t) dt 


E 1 
M CON EE IR (13) 


u 
h f 
dt 
ti 


where t; and t? are travel times from the dissemination site to particu- 
lar arcs. From equations (11) and (8), 


t2 » 
f [hi +К12 (te - t9] dt i" 1 aoe !) 
t a = 
1 А 2 1 
Gp ge з= ы ысышы шл cer (14) 


t 
f 2 at (n +1) (h2 - hj) 
ti 


where, again, the subscripts indicate particular arcs. 


3 ST. LOUIS DISPERSION STUDY II 


For particular arcs 


x, - 5 

= 2 

a ve 
2 1 


where x is the mean distance from dissemination site to an arc and t 
is the corresponding travel time based on the appropriate value of 9). 
Combining equations (14) and (15), 


2 
бый Ош жы arn 
922 T 


For each evening experiment, the coefficients a and n were first 
determined from mean wind profiles measured at the downtown TV 
tower or by pilot balloon ascents, after which the h for the appropriate 
sampling arcs was computed from equation (12). А value for ўр cor- 
responding to each h was calculated from equation (11), and travel 
times between arcs were determined from equation (16). 


TABULATION OF RESULTS 


Symbols used in the foregoing calculations are listed in Table 1, 
and dispersion data resulting from the calculations are summarized 
in Table 2. In addition, in Table 2 values of an "equivalent! Gaus- 
sian ez, i. e., 2Һ/,/ От, are presented for evening experiments іп 
order to permit a single analysis of dispersion in the vertical. 


LOSS OF TRACER MATERIAL 


Since disseminations were from near ground level, the loss of 
tracer material due to deposition or fallout was thought to be a possi- 
bility. For the analyses reported herein, such losses would most 
noticeably affect estimates of vertical dispersion parameters and of 
axial concentration. A special study (Vaughan and McMullen, 1968) 
was conducted to assess the magnitude of deposition or fallout that 
occurred during the St. Louis tracer experiments under various mete- 
orological and ambient air pollution conditions. The assessment was 
based on statistical comparisons of percentages of fluorescent parti- 
cles in class-size groupings between samples collected at the various 
downwind distances and between these collected samples and the stan- 
dard "lot" disseminated. The findings indicated that in general signifi- 
cant losses of tracer material did not occur in transit. An apparent 
loss of small particles in the size range less than 1. 5 microns did 
occur with distance under stable meteorological conditions. This loss 
seemingly resulted from background obscuration of collected particles 
during prolonged sampling times at greater distances, and not from 
deposition or fallout. In the St. Louis study, particles with less than 
1. 5-micron diameters constituted from 20 to 30 percent of the total 
mass released during each experiment. 
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In addition, tracer material collected by samplers may not be 
measured because of loss of fluorescence, particularly when exposed 
to high humidity and ultraviolet radiation. A recent study (Grinnel, 
1965) has demonstrated that the loss in fluorescence of the type of 
tracer material utilized in St. Louis is no more than 5 percent. 


EFFECTIVE TRANSPORT WIND SPEED 


As previously discussed, for daytime experiments a single effec- 
tive transport wind seemed to exist and was represented best by the 
mean tetroon trajectory, when the tetroon was found in the airflow 
dispersing the tracer cloud. It was assumed that during the more 
stable conditions of the evening experiments the vertical profiles of 
wind and tracer distribution were interrelated, resulting in a transport 
wind increasing with vertical spreading of the tracer cloud, and thus 
with travel time. 


Instrumental measures of wind speed over tracer dissemination 
periods were qualitatively evaluated in terms of mean transport wind 
speeds. Values computed from average pilot balloon observations were 
of limited use, since deviations from the average (at a particular level) 
usually were a significant fraction of that average. Wind speed meas- 
urements at the tracer dissemination sites provided insufficient data 
to permit meaningful comparisons with other wind data. Wind speeds 
measured at about 60 feet above ground at the peripheral sites and at 
the 127-foot level of the downtown tower were consistently lower than 
both the daytime and evening effective transport wind speeds. Meas- 
urements at the 255-foot level of the tower, during the limited number 
of experiments for which this wind equipment was operating, usually 
were of the appropriate magnitude for both daytime and evening effec- 
tive transport winds. Although measurements at the 459-foot level of 
the tower usually provided better estimates of daytime transport wind 
speeds than did measurements at the 255-foot level, they generally 
overestimated evening transport wind speeds. 


Table 1. LIST OF SYMBOLS DESCRIBING DISPERSION DATA 


x Downwind distance from dissemination site 

t Travel time from dissemination site 
су Cross-wind standard deviation of tracer material 
Cz Effective vertical standard deviation of tracer 

material 
h Height of uniformly mixed layer 
(x/Q)p Relative axial concentration 

а Effective transport wind speed 

u Wind speed at height z 

2, Height 

a,n Empirical constants 
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METEOROLOGICAL INDICES OF TURBULENCE 


The calculated values of dispersion parameters were related to 
several readily measured or derived meteorological indices of turbu- 
lence. These indices included stability classes, gustiness classes 
based on patterns of wind direction fluctuations, and a classification 
based jointly on measures of wind direction fluctuations and conditions 
of vertical stability. Values of the indices were determined for semi- 
rural or peripheral locations as well as for urban locations. This ap- 
proach allows at least qualitative inferences concerning the applicability 
of the various indices and of station locales in describing dispersion 
over urban areas. 


PASQUILL-TURNER STABILITY CLASSES 


In the scheme devised by Pasquill (1961) and later slightly modi- 
fied by Turner (1964), the diffusive ability of the lower atmosphere is 
represented by stability classes determined jointly from estimates of 
net radiation and wind speed . Determinations of these stability classes 
over periods of tracer dissemination were made for three wind meas- 
urement locations in the St. Louis area: the lower level (127 feet above 
ground) on the TV tower, the peripheral site west of downtown, and the 
Weather Bureau Airport Station at Lambert Field in comparatively rural 
surroundings. The net radiation estimates were based on solar altitude, 
obtained from the Smithsonian Meteorological Tables (1963), and cloud 
cover and ceiling height, taken from records of hourly observations at 
the airport. 


The use of Pasquill-Turner stability classes for peripheral and 
urban locations may, of course, be questioned, since the scheme was 
devised primarily for application over relatively open country. Also, 
information on cloud cover at the airport may not always be appropri- 
ate for the other locations. 


MODIFIED “BROOKHAVEN” GUSTINESS CLASSES 


In the gustiness classification scheme due to Singer and Smith 
(1953) at Brookhaven National Laboratory, the diffusive ability of the 
lower atmosphere is described by the range and rapidity of horizontal 
wind direction fluctuations; a less detailed scheme had been designed 
previously by Giblett et al. (1932) for use at Cardington. The range 
furnishes an estimate of the turbulent velocity component in the cross- 
wind direction and hence of cross-wind dispersion. The rapidity of 
the fluctuations also indicates the general degree of vertical stability 
and thus of dispersion in the vertical. 


The basic Brookhaven gustiness classes were used in St. Louis, 
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but different class limits on the ranges of wind direction fluctuations 
were adopted. Although the St. Louis wind sensors were the same type 
as those used at Brookhaven, they were located near larger roughness 
elements and usually nearer the ground. In addition to the urban (lower 
level on TV tower) and peripheral sites discussed earlier, wind data 
for the upper level (459 feet above the ground) of the TV tower were 
utilized also. With these additional data, comparisons can be made 
between wind direction fluctuations at 127 feet, slightly above nearby 
building tops, and those at 459 feet, well above the buildings. For each 
site, range limits for the specific classes were derived from appropri- 
ate analog traces obtained from a climatological analysis of the (analog) 
records and were, without exception, greater than those found by Singer 
and Smith (1953) at the Brookhaven site. 


HORIZONTAL WIND DIRECTION FLUCTUATIONS 
AND CONDITIONS OF VERTICAL STABILITY 


This joint meteorological index employs the standard deviation of 
horizontal wind direction fluctuations, 90, and a direct estimate of 
vertical stability expressed as a gradient Richardson number. Indices 
based partially or totally on these elements have been applied by such 
investigators as Cramer (1959), Fuquay et al. (1964), and Slade (1965) 
to organize dispersion data from other experimental programs. 


Estimates of % over the tracer dissemination periods were made 
from frequency histograms of azimuth angle derived from records of 
the analog wind traces. Since the chart speed of wind recorders was a 
slow 3 inches per hour, portions of the analog traces were occasionally 
"painted" (1.е., the more commonly occurring azimuth directions). 
Frequencies for the ''painted'' azimuth angles were extrapolated from 
those of the "non-painted'' angles under the assumption that the fre- 
quency distributions of azimuth angle were Gaussian or resulted from 
sums of Gaussian distributions. 


Comparisons of these values of % with the extreme, third highest, 
and fifth highest ranges of wind direction fluctuations over the tracer 
dissemination periods showed statistically linear relationships for each 
range for each meteorological site. Linear correlation coefficients 
varied from 0.8 to slightly over 0.9, the higher correlations generally 
occurring for the peripheral and upper-level urban locations. Some- 
what better correlations were attained with the third highest range than 
with the extreme or fifth highest ranges. Average values of the ratio 
of range tovg for the extreme, third highest, and fifth highest range 
were, respectively, about 7.5, 6.0, and 5.5. The value of 7.5 for the 
extreme range is slightly different from the 6. 0 found by Markee (1963) 
and Slade (1965) for other locations and exposures. 


Estimates of stability were based on wind and temperature meas- 
urements at the downtown TV tower or temperature measurements from 
the downtown radiosonde ascents. Lapse rates, stability ratios, and 
various forms of the gradient Richardson number for several layers in 
the lower atmosphere were considered. Thorough evaluation of these 
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parameters was difficult because of incomplete data. А limited ap- 
praisal, however, indicated that the most meaningful parameter was 
the "bulk" Richardson number Rip. Fortunately, adequate data were 
available to compute this value for most experiments. The Rip was 
based on temperature and wind measurements at the 127- and 459-foot 
levels of the TV tower. Temperature information from the radiosonde 
ascents was substituted when the corresponding tower data were miss- 
ing. Following Lettau (1957), the Rip may be defined: 


fat/az+ Тај 22 


1 = _Ё- 
"ig. wp 72 
V 
where: 
g = acceleration of gravity 
AT - temperature difference between top and bottom of the 
layer 


AZ = height difference between top and bottom of the layer 
T = mean absolute temperature through the layer 

Га = dry adiabatic temperature lapse rate 

Z = height of upper anemometer 

V = mean speed for anemometer at height 2, 


Originally, a measure of a transport wind speed, ū, was also in- 
corporated into the joint index; its utilization as an additional parameter 
or in the form 9, u generally increased data scatter. It should be noted, 
however, that the tracer experiments were not usually conducted under 
very strong or very light wind speed conditions, resulting in a reduc- 
tion in the overall and perhaps within-class range of 1. Dispersion in 
the cross-wind and in the vertical were better ordered by og and Rip 
jointly than by either alone. Axial concentrations also were better re- 
presented by this joint classification than by og or Rip alone. 


TABULATION OF RESULTS 


For each of the experiments, values of the meteorological indices 
at the locations specified in the preceding sections are presented in 
Table 3. "Values of Rip denoted with superscripts were those that were 
inconsistent with all other coincident measures of stability, e.g., Pas- 
quill-Turner classes or gustiness classes. Without exception, these 
inconsistencies were for cases in which the radiosonde temperature 
data were substituted for missing TV tower data. 
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RELATION OF DISPERSION PARAMETERS 
TO METEOROLOGICAL INDICES OF TURBULENCE 


Values of ту» ©, (Х70), and (x/Q)pu were plotted against down- 
wind distance, x, and travel time, t. For daytime experiments, 8 was 
based on tetroon speeds; for evening experiments, it was Ti, (see Verti- 
cal Parameters). The dispersion data points were grouped according 
to values of the meteorological indices determined for each of the pre- 
viously mentioned locations. Best-fit curves through the data groupings 
were approximated by eye in view of the small sample size. Linear 
curves were used since, for the arc distances and travel times repre- 
sented in the data, few indications of discernible non-linearity were 
noted. Data points іп correlograms involving e, and axial concentra- 
tion were ignored in the determination of best-fit lines when a re- 
striction to vertical mixing (e.g., a temperature inversion aloft) was 
thought to be affecting the concentration. For such cases the appro- 
priate data in Table 3 are denoted with а superscript a. Information 
concerning the depth of a mixing layer was obtained from the downtown 
radiosonde ascents and less frequently from the regular Weather Bu- 
reau rawinsonde ascents at Columbia, Missouri, or Peoria, Illinois. 
The determination that a particular mixing depth affected the dispersion 
data was made from a consideration of the mixing depth, the growth 
rate of c, and the travel time or downwind distance. Unfortunately, 
sufficient definitive information was not available to permit a detailed 
analysis of these cases. 


Correlograms of оу, 77, (X/Q)p, and (X/Q)pu versus downwind 
distance in terms of øg and Rig for the lower-level urban location are 
shown in Figures l through 4. Similar plots in terms of travel time 
are presented as Figures 5 through 8. Also, for this urban site, cor- 
relograms of су, 92, апа (Х/О)р versus downwind distance are given 
in Figures 9 through 11 as functions of Pasquill-Turner stability classes; 
as functions of modified Brookhaven gustiness classes they are given 
in Figures 12 through 14. Relations (not shown) in terms of travel time 
for these other indices were similar to those in Figures 5 through 8. 


The scatter of dispersion data points about best-fit lines was 
greater in terms of meteorological indices computed at peripheral or 
semirural locations than for the same indices at the lower-level urban 
locations. Relationships for the dispersion data in terms of indices at 
the upper-level urban location were similar to those in terms of the in- 
dices for the corresponding lower-level location. 


In Figures 1 through 14 the following additional features are wor- 
thy of specific note: 
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1. The ordering of data was poorest for the close-in sampling 
arcs at which multi-peaked or otherwise complex distributions in the 
cross-wind occurred in greater proportion than at more distant arcs. 
As previously noted, channeling of the airflow due to various locally 
induced circulations may be largely responsible for the appearance of 
such distributions. 


2. The scatter of data points about best-fit lines was much 
greater in о, than in the other correlograms. Since the actual vertical 
distributions of tracer material were largely unknown, the relative 
magnitude of such scatter was not altogether unexpected. 


3. When plotted according to travel time, data were grouped best 
into only ''stable"' and "non-stable'' regimes, except for о„, which was 
ordered about as well by travel time as by downwind distance. 


4. Inan earlier preliminary analysis of the dispersion data, 
Pooler (1966) found that values of 9.,, ordered only according to day- 
time and evening categories, were grouped about equally well in terms 
of downwind distance and of travel time, as shown in Figures 1 and 5. 
In these figures daytime is indicated with a very few exceptions by 
those points for which Rip < 0.01. As illustrated in Figures 2 and 6, 
the ordering of 9, in terms: of daytime and evening categories is better 
in terms of travel time than of downwind distance. 


5. Values of c, and axial concentration that were considered to 
be significantly affected by a restrictive mixing layer were usually 
smaller and larger, respectively, than their non-affected counterparts 


under similar meteorological conditions. 


6. Onlya slight improvement in ordering of data points was 
achieved through normalizing (Х/0)р by u. Since, however, as was noted 
earlier, the range of likely ü's is not fully represented in the data here, 
the generality of this result and of that involving the ordering of v, with 
distance and time may require further verification. 


7. The scatter of data points about best-fit lines was considerably 
greater for Pasquill-Turner stability classes than for the other indices. 
For the modified Brookhaven gustiness classes and the joint classes of 
сө and Rip, the scatter of data points about best-fit lines was quite 
similar. 
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COMPARISONS WITH OTHER EXPERIMENTAL 
PROGRAMS 


The results (i.e., the derived dispersion relationships) of the 
St. Louis tracer experiments сап be compared with those of past ''open 
country'' programs for which the dispersion parameters are described 
by comparable meteorological indices. Specifically, the St. Louis 
results based on values of ¢gand Rip (Figures l and 3) are related to 
those of the Green Glow - Series 30 experiments (Fuquay et al. , 1964); 
the results according to Pasquill-Turner stability classes (Figures 9 
and 10) are compared with those of several programs as summarized 
by Pasquill (1961), modified slightly by Gifford (1961) and presented by 
Hilsmeier and Gifford (1962); and the St. Louis results according to 
modified Brookhaven gustiness classes (Figures 12 and 13) are related 
to those of the Brookhaven experimental program (Singer and Smith, 
1966). These comparisons (Green Glow, Pasquill-Turner stability, 
and Brookhaven gustiness) are illustrated in Figures 15 through 17. 


The analysis by Fuquay et al. (1964) is in terms of travel time. 
In view of the St. Louis experimental results, the Green Glow - Series 
30 data were analyzed in terms of downwind distance (Figure 15). Of 
the meteorological parameters considered (сё, ч, cg 4, and a gradient 
Richardson number), dispersion in the cross wind for Green Glow - 
Series 30 was best described by о, alone. Peak exposure (axial con- 
centration) was best represented by joint measures of ¢ and the Rich- 
ardson number; peak exposure normalized by a ü was not ordered as 
well. Data reported for the 12. 8- and 25. 6-kilometer arcs were con- 
sidered insufficient to facilitate a detailed analysis with respect to the 
meteorological parameters and hence were not used here. 


In Figure 16 the results of the St. Louis experiments are compar- 
ed with results of similar work in other areas. The grouping of the 
St. Louis, Johnstown, and Fort Wayne data herein was by Pasquill- 
Turner stability classes (Turner, 1964), which are merely more ob- 
jective expressions of the Pasquill classes in terms of readily avail- 
able meteorological variables. The solid curves are based upon data 
ordered by the Pasquill stability classes (Pasquill, 1961); those in 16a, 
and 16b. were originally presented by Gifford (1961) and in 16с. by 
Hilsmeier and Gifford (1962). 


Particularly in Figures 16a and 165, note that the best-fit lines 
for су and о; for the St. Louis data are everywhere larger than those 
for the open-country counterparts. When extrapolated to longer down- 
wind distances, the best-fit lines for Ty approach their counterparts; 
those for % generally approach counterparts of stability one class 
higher. In Figure lóc, the dilution of (x/Q)gu in St. Louis over that in 
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Open country is greatest near the source and is greater for stable than 
for unstable meteorological conditions, as expected (Figures l6a and 
165). Overall relations similar to those in Figure 16 are generally 
noted in Figures 15 and 17 despite the differences in terrain, sensor 
exposures, and response characteristics of equipment used in the 
various studies. Some of the apparent differences in Figure 17 can 
also be accounted for by the fact that Pasquill-Turner class A's and 
F's, which constitute a sizable fraction of the Brookhaven Вә and D 
class members, respectively, were seldom represented in the data 
collected in the St. Louis experimental series. 


Values of dispersion parameters obtained in fluorescent particle 
tracer experiments in Johnstown, Pennsylvania (Smith, 1967), and 
Fort Wayne, Indiana (Csanady et al., 1967), are also presented in 
Figure 16. The values are similar to those obtained in St. Louis 
under the same overall meteorological conditions. In Johnstown, dis- 
seminations were from a low-level point source in the urban area; in 
Fort Wayne they were from a 90-meter-high, cross-wind line source 
located 1 mile upwind of the urban area. Values оѓо, for this elevated 
line source that were reported for the 1-mile downwind distance (i.e., 
at the upwind edge of the urban area) are not shown here; a direct com- 
parison of the other values, particularly for shorter distances, with 
those of low-level sources can be questioned since the effects of local 
obstructions may not have become significant until the tracer plume 
reached near ground level. 


For the St. Louis experiments a ratio of peak to mean concentra- 
tions (P/M) was computed for each case in which the appropriate se- 
quential sampler was on or quite near the mean centerline of a tracer 
cloud. The ratios were formed from the peak and arithmetic average 
(mean) values of concentration of the time-concentration histograms; 
intervals of beginning and ending 'dribble'" were excluded. Plots of 
P/M as a function of downwind distance and travel time showed little 
variability, even when data were separated according to stability. 
Gifford (1960) presents similar results for ground-level emission 
sources. 


A plot of P/M as a function of the ratio of averaging time, t,, to 
sampling time, t,, is shown in Figure 18. Averaging time is the time 
interval required for passage of the tracer cloud, excluding beginning 
and ending dribble. Sampling time is the finite time increment, 
usually 1 or 2 minutes, over which the peak concentration was meas- 
ured. In this figure, P/M values show a very slight tendency to in- 
crease with increasing t,/t,. The values of P/M vary from about 2 to 
6. Gifford (1960) shows similar values for ground-level sources; 
Singer et al. (1963) found that values decreased in proportion to in- 
creasing density of vegetation near receptors. 
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INITIAL DIMENSIONS OF TRACER CLOUD 


The effects of a building upon a plume generated in its vicinity 
have been studied extensively in wind tunnel experiments (e.g., Halit- 
sky, 1963). More recently, a full-scale investigation was conducted 
. in the atmosphere (Dickson et al., 1967). Essentially, these studies 
show that the primary effect is an enhancement close to the source of 
the size of the plume, hereafter referred to as the initial size of the 
plume. 


Pasquill-Gifford (Gifford, 1961) "С" stability curves (e.g., see 
Figure 16) for о, and о, and revised curves allowing for virtual travel 
distances corresponding to the specified initial cloud dimensions are 
presented in Figure 19. То facilitate comparisons, the best-fit "С! 
lines for the St. Louis tracer data (Figure 16) are also shown. The 
data in Figures 16 and 19 along with those in 15 and 17 indicate that the 
effect of the city on dispersion in the cross-wind is similar to the re- 
cognized effect of a single building near the source, although more 
complex. When the cloud becomes much larger than eddies created by 
the local obstructions, the extent of the dispersion, or rather 9, close- 
ly converges to that associated with open-country flow; the initial en- 
hancement of plume size is presumably dependent on the width and 
height of the obstructions. A large initial dimension of the cloud and 
considerably enhanced dispersion is indicated in the vertical. The 
enhancement in the rate of dispersion over that in open country is 
greater in stable than in unstable meteorological conditions and is pre- 
sumably due largely to enhanced convective activity over the urban 
environs. 

In Figure 19, estimates of 50 to 60 meters for an initial о, and 
20 to 30 meters for an initial c, for the somewhat diverse St. Louis 
locations appear reasonable. Relationships for the other stability 
classes are similar to those for "С." A Gaussian plume from a 
ground-level point source with an initial depth (2. 15 «,, where plume 
edge - 1/10 (x/Q)p) equivalent to the average building height of 40 to 
50 meters in downtown St. Louis would have а с, value of approxi- 
mately 20 meters. For a similar plume with an initial width (4.3 9.) 
equivalent to the length of a typical city block (160 meters), ao ot 
roughly 40 meters would result; if the plume were oriented with the 
cross -wind direction diagonally across such a block, the initial value 
of ø would increase to about 60 meters. Particularly in the down- 
town and older residential areas of cities, streets usually constitute 
the only significant relief between structures. 
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Cross -wind standard deviation of tracer material as a function 
of downwind distance in terms of standard deviation of wind 
direction fluctuations (26) and bulk Richardson number (Rig). 
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Figure 2. Effective vertical standard deviation of tracer material as a 
function of downwind distance in terms of standard deviation 
of wind direction fluctuations (c9) and bulk Richardson 
number (Rig). 
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Figure 3. Relative axial concentration as a function of downwind distance 
in terms of standard deviation of wind direction fluctuations (96) 
and bulk Richardson number (Rig). 
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Figure 4. Normalized relative axial concentration as a function of down - 
wind distance in terms of standard deviation of wind direction 
fluctuations (7g) and bulk Richardson number (Rig). 
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Figure 5. Cross - wind standard deviation of tracer material as a function 
of travel time in terms of standard deviation of wind direction 
fluctuations (сө) and bulk Richardson number (Rig). 
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Figure 6. Effective vertical standard deviation of tracer material as a 
function of travel time in terms of standard deviation of wind 
direction fluctuations (95) and bulk Richardson number (Rig). 
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Figure 7. Relative axial concentration as а function of travel time in 


terms of standard deviation of wind direction fluctuations (00) 
and bulk Richardson number (Rig). 
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Figure 8. Normalized relative axial concentration as a function of travel 
time in terms of standard deviation of wind direction fluctuations 


(25) and bulk Richardson number (Rig). 
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Figure 9. Cross -wind standard deviation of tracer material as a function 
of downwind distance in terms of Pasquill - Turner stability 
classes. 
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Figure 10. Effective vertical standard deviation of tracer material as a 
function of downwind distance in terms of Раѕдиі!! - Turner 
stability classes. 
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Figure 11. Relative axial concentration as a function of downwind dis- 
tance in terms of Разаш! - Turner stability classes. 
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Figure 12. Cross -wind standard deviation of tracer material as a 
function of downwind distance in terms of modified Brook- 
haven gustiness classes. 
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Figure 13. Effective vertical standard deviation of tracer material as 
a function of downwind distance in terms of modified Brook- 
haven gustiness classes. 
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Figure 14. Relative axial concentration as a function of downwind distance 
in terms of modified Brookhaven gustiness classes. 
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Figure 16 (а and- b). Comparison of the results of the St. Louis and other urban tracer 
experiments with those summarized by Pasquill and Gifford (1961). 


Initial Dimensions of Tracer Cloud 39 


Note: 

The grouping of the St. Louis, Johnstown, and Fort Wayne data in Fig- 
ure 16 was by Pasquill-Turner stability classes (Turner, 1964), which 
are really only more objective expressions of the Pasquill stability 
classes (Pasquill, 1961) in terms of readily available meteorological 
variables. The solid curves are based upon data ordered by the Pas- 
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Figure 16 (с). Comparison of the results of the St. Louis tracer experiments with those 
summarized by Pasquill and Gifford (1961). 
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Figure 17. Comparison of results of St. Louis tracer experiments with those 
of Brookhaven experiments (Singer and Smith, 1966). 
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Figure 18. Ratio of peak to mean concentration as a function of ratio 
of averaging time to sampling time. 
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Initial Dimensions of Tracer Cloud 


CONCLUSIONS AND RECOMMENDATIONS 


Results of the St. Louis Dispersion Study reported here support 
the following conclusions regarding dispersion of airborne material 
emitted from low-level point sources in urban areas: 


1. Dispersion can be readily described by commonly utilized 
meteorological indices of turbulence. The more detailed indices for 
urban locations appeared to be most representative. 


2. In terms of the meteorological indices of turbulence, cross- 
wind dispersion is better described as a function of downwind distance 
than of travel time, whereas vertical dispersion is described about as 
well by travel time as by downwind distance. It should be noted, how- 
ever, that Pooler (1966) found, for an overall classification of experi- 
ments simply as daytime or evening, that cross-wind dispersion was 
expressed about as well in terms of travel time as of downwind dis- 
tance, whereas vertical dispersion was expressed better in terms of 
travel time. 


3. The urban area affects cross-wind dispersion primarily by 
enhancing the size of the initial plume. When the plume becomes much 
larger than the size of eddies created by the local obstructions, the 
extent of the dispersion approaches that associated with flow over open 
country. In the vertical, significantly enhanced dispersion as well as 
a large initial spread of the plume result; the enhancement in the rate 
of dispersion over that in open country is somewhat greater for stable 
than for unstable meteorological conditions and, presumably, is due 
largely to enhanced convective activity over the urban environs. 


4. Restrictive mixing layers (e.g., inversion aloft) can signifi- 
cantly alter the values of affected vertical dispersion parameters and 
concentrations of airborne material near the surface. 


Dispersion from low-level sources in urban areas for downwind 
distances of less than about 1/2 mile is conjectural. Here, the effects 
of local roughness elements should be most pronounced. Specifically, 
initial cloud dimensions as affected by factors such as building width 
and height and building density should be ascertained, and the effects 
of complexes rather than single structures should be catalogued. 


The effects of urban areas on dispersion of plumes from elevated 
sources, particularly at the shorter downwind distances, may differ 
markedly from those suggested for lower-level sources, and hence may 
require independent investigation. As mentioned earlier, the effects 
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of local obstructions may not become significant until a plume nearly 
reaches the ground. At this point the plume may also be so large that 
the locally created eddies are not significantly effective in dispersing 
it. In addition, as the effective source height increases, different 
stability regimes about which little is known are encountered. 


The analyses reported herein concerning data collected during 
the St. Louis Dispersion Study are not considered to be exhaustive. It 
is hoped that the body of data can serve a useful purpose for the appli- 
cation of more advanced analytical techniques, which are certain to be 
developed by meteorological science in the future. 
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